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Abstract 

We consider the problem of reconstruction of planar domains from their moments. 
Specifically, we consider domains with boundary which can be represented by a 
union of a finite number of pieces whose graphs are solutions of a linear differen- 
tial equation with polynomial coefficients. This includes domains with piecewise- 
algebraic and, in particular, piecewise-polynomial boundaries. Our approach is 
based on one-dimensional reconstruction method of [5] and a kind of "separation 
of variables" which reduces the planar problem to two one-dimensional problems, 
one of them parametric. Several explicit examples of reconstruction are given. 

Another main topic of the paper concerns "invisible sets" for various types 
of incomplete moment measurements. We suggest a certain point of view which 
stresses remarkable similarity between several apparently unrelated problems. In 
particular, we discuss zero quadrature domains (invisible for harmonic polyno- 
mials), invisibility for powers of a given polynomial, and invisibility for complex 
moments (Wermer's theorem and further developments). The common property 
we would like to stress is a "rigidity" and symmetry of the invisible objects. 
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1 Introduction 



In this paper we continue our study of nonlinear problems of reconstruction 
of multidimensional objects from incomplete collection of integral measure- 
ments. The paper has two parts, closely related but different in their goals. 
In the first part we present a method of reconstruction of planar domains 
of a certain special class from finite collections of moments. In the second 
part we discuss the structure of sets and functions "invisible" for a certain 
collection of moment measurements. 

In more detail, the object we would like to reconstruct is a 2-dimensional 
finite domain GcK 2 which we assume to belong to a certain finite-dimen- 
sional family G\ specified by a finite number of discrete and continuous 
parameters A. Specifically, we shall assume that the boundary of G is a 
union of a finite number of pieces whose graphs are solutions of a linear 
differential equation with polynomial coefficients. 

The measurements are represented by finite collections of the moments 
m Q>( 3 of the characteristic function xg(x,v) of the domain G: 



Our main problem is to provide an explicit (and potentially efficient) recon- 
struction method, and, in particular, to estimate a minimal possible set of 
these moments sufficient for unique reconstruction of the domain G. 

Similar inverse problems have been intensively studied, including recon- 
struction from their moments of polygons, of quadrature domains, of certain 
"dynamic" semi-algebraic sets, see [7J El [T7J |20j [33] and references therein. 
In a more general context the problem of domain reconstruction from its mo- 
ments appears as a part of broad field of inverse problems in Potential Theory 
(see, for example, [HE])- Rather similar questions arise in reconstruction from 
tomography measurements ([HI [191 EI])- 

Our approach is based on one- dimensional reconstruction method of [5] 
applicable to piecewise continuous functions satisfying on each continuity 
interval a linear differential equation with polynomial coefficients. Then we 
use a kind of "separation of variables" which reduces the planar problem to 
two one-dimensional problems, one of them parametric. 

We expect that a reconstruction method for piecewise-smooth functions 
given in [9] can be extended in a similar way also to planar and higher 
dimensional piecewise-smooth functions. 
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The second part of the present paper is devoted to "invisible sets" for var- 
ious types of incomplete measurements. Here we do not provide new results 
(besides a couple of examples), but rather suggest a certain point of view 
which stresses remarkable similarity between several apparently unrelated 
"moment vanishing" problems. In particular, we discuss zero quadrature 
domains (invisible for harmonic polynomials), invisibility for polynomials 
annihilating other partial differential operators, invisibility for powers of a 
given polynomial, and invisibility for complex moments (Wermer's theorem 
and further developments). In all these cases we stress a common property 
of "rigidity" and symmetry of the invisible objects. 

2 One-dimensional case 

Reconstruction problem in dimension one has been settled in a pretty sat- 
isfactory way for many important finite-dimensional families of functions. 
This includes linear combinations of shifts of known functions, signals with 
"finite rate of innovation", piecewise D-finite functions which we use below, 
piecewise-smooth functions, and many other cases (see [5], HI [HI [35] and 
references therein). 

2.1 Piecewise D-finite reconstruction 

Let g(x) be a function with a support [a, b] C E 1 , satisfying the following 
condition: there exists a finite set of /C + 2 points a = £o < £i < • • • < 6c+i = 
b, such that on each segment [£n>£n+i]j n — 0, 1, . . . 1C, the function g(x) is 
continuous and satisfies there a linear differential equation 




with polynomial coefficients p n> j(x) = Yl a n,i,j%\ Pnj ^ on [a, b}. At the 

i=0 

points £ n the function g(x) may have jumps. Such functions are described 
by a finite collection of discrete and continuous parameters, and they are 
called piecewise D-finite. Without loss of generality (at least theoretically) 
we can assume that all the operators D n = D are the same. In particular, 
piecewise-algebraic, and, specifically, piecewise-polynomial functions belong 



2 



to this class. In the last case the differential operator is D 



_ d N + 1 

~ dx N + 1 ' 



where 



N is the maximal degree of the polynomial pieces of g. 

It was shown in j5] that the collection of "discrete" parameters /C, N, 
{k n j}, together with a sufficiently large collection of the moments 



determine uniquely any D-finite function g(x) with all the points £o> • • • C/c+i 
of its possible discontinuity, as well as the coefficients of the differential op- 
erator D. 

For piecewise-polynomial functions the number \x of the moments required 
for reconstruction, depends only on the discrete data: the number of jumps 
K and on the maximal degree N of the pieces. It is shown in [5] that in 
piecewise-polynomial case 



A similar, but more complicated expression for /x can be written in 
piecewise-algebraic case. However, for general D-finite functions, with re- 
spect to a general second (and higher) order differential operators D, the 
number \i may depend also on specific coefficients of D. 

Let us give a very simple example of this latter phenomenon. Let L n (x) be 
the n-th Legendre polynomial, defined 

2^~~t T™\( x2 ~ -Q n ]- Legendre 
polynomials are pairwise orthogonal on [—1, 1] and they satisfy the second 
order Legendre differential equation — x 2 )-^L n (x)] +n(n + l)L n (x) = 0. 
Since Lj(x), j < n, form a basis of the space of all polynomials of degree n we 
conclude that L n (x) is orthogonal to 1 . Hence the moments 

rrij(L n ) = J_ 1 xi L n (x)dx vanish for j = 0, 1, . . . , n — 1. We conclude that 
a D- finite function L n (x) on [—1, 1] cannot be reconstructed from less than 
n + 1 its moments. So \x above depends not only on the order and degree of 
the Legendre operator, but also on a specific value of the parameter n in it. 

Notice that the leading coefficient of the Legendre equation vanishes at 
both the endpoints —1,1 of the interval. 

The reconstruction procedure described in [5] consists of solving certain 
linear and non-linear algebraic equations whose coefficients are expressed 
through the moments m a . These equations have a very specific structure 




g(x) ■ x a dx, a = 0, 1, 2, . . . , \i 



/i = /i(/C, N) = max <^ 2(iV + l)K - 2, (/C + 1){N + 1) 
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which we illustrate in the next section with the simplest example of the 
classical "Prony system". The last step requires also finding a basis of the 
solution space of the differential equation Dg = 0. 

We shall apply below this reconstruction procedure, referring to it as to 
Procedure 1. 

2.2 Prony system 

Prony system appears as we try to solve a very simple version of the shifts 
reconstruction problem. Consider F(x) = Y^j=i a j3( x ~ x i)- We use as 
measurements the polynomial moments 



This infinite set of equations is called Prony system. It can be traced at least 
to R. de Prony (1795, [32]) and it is used in a wide variety of theoretical 
and applied fields. See, for example, [6J [35] and references therein for a very 
partial list, as well as for a sketch of one of the solution methods. This method 
requires 2N equations from (2.2). It allows first to find the number of nonzero 
coefficients aj. Then dj and Xj are found via solving first a Hankel-type linear 
system of equations with coefficients formed by the moments mi, interpreting 
the solution as the coefficients of a certain polynomial, and finding all the 
roots of this polynomial. 

We shall apply below this Prony solution procedure in our specific situa- 
tion, referring to it as to Procedure 2. 

3 Main result 

We assume that the domain G C M 2 to be reconstructed has a "D-finite 
boundary" . More accurately, we have the following definition: 




After substituting F into this integral we get m n = J Ylj=i diS(x — Xj)x n dx 
Yuf=i a j x ]- Considering Oj and Xi as unknowns, we obtain equations 



N 




(2.2) 
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• • • 



6c &+i = b 



Figure 1: Schematic representation of the domain G. 



Definition 3.1 A compact domain G CR 2 is called D-finite if the boundary 
dG is a union of k segments Sj with the following property: there exists a 
linear differential operator D of the form (2.1) and with the leading coefficient 
not vanishing for x in the projection of G, such that each Sj is the graph of 
a function y = ipj(x) satisfying Dipj = 0. 



In particular, if each Sj is a graph of an algebraic function y = ipj (x) and all 
the branches of these algebraic functions are regular over the projection of 
G, then G is D-finite. The simplest but still important example of D-finite 
domain is when all tpj are polynomials. 

We do not restrict the topological type of G — it may have "holes" . 

Before we formulate the main result, let us introduce some notations. Let 
[a, b] be the projection of G onto the x-axis, and let a = £o < £i < • • • < 
6c+i = b be all the projections of the endpoints of the segments Sj of the 
boundary dG (see Figure 1). Certainly, K, < k, while the maximal number 
of the intersection points of dG with vertical lines is at most k — 1. 
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Write the linear differential operator D in Definition 3.1 as 

Dg{x)=Y^p i {x)^{x) i (3.1) 



3=0 



where Pj(x) = a ijx\ 



i=0 



Theorem 3.1 Any D-finite domain with the discrete parameters k, N, kj can 
be uniquely reconstructed from a collection of double moments 

m a ,/3 = J J 2 Xg(x, y) ■ x a y^dxdy, 0<a< M(k, N, kj), < (3 < 2(k - 1). 

(3.2) 

The reconstruction procedure requires solving of certain linear and non-linear 
algebraic equations whose coefficients are expressed through the moments 
m a fi, and solving equation Du = with specific numerical coefficients found 
in previous stages. 

Proof: Denote by Aj the interval j — 0, . . . ,/C. Over each Aj the 

domain G is a union of Sj < |(« — 1) strips (f> < y < 0^, , I = 1, . . . , Sj, 

(see Figure 1). We have 



«5 K, p 

m«,/3 = / x a Vp(x)dx = J2 x a ^^(x)dx, (3.3) 

J a j_Q J A, 



where for x G Aj 



*f>A*) = r ^ v"xg(x, y)dy = -L- f^[$\x) - <ffi\x)]. (3.4) 

The first conclusion is that for each f3 > the function typ is piecewise D- 
finite. Indeed, on each interval Aj, j = 0, ...,/C, ^ p = \P pj is a linear 
combination of the /3-th powers of the functions 0,0^, Z = 1, . . . , Sj which, 
by the assumption on G, satisfy Dip = 0. Hence ^fp (i.e. each satisfies 
another linear differential equation with polynomial coefficients Dptyp = 0. 
The operator Dp depends only on D, in particular, its order and degree 
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depend only on the order and degree of D, and it has no singularities on 
[a, b], if D possesses this property. 

We can find one-dimensional moments of ^ p via (3.3): 



Now we apply one-dimensional Procedure 1 from Section |2TT| and reconstruct 
typ from the moments m a ,p, a = 0,1, . . . , [ip < M(k, N, kj) for each (3 = 
0, 1, . . . , 2(k — 1). The reconstruction procedure requires solving of certain 
linear and non-linear algebraic equations whose coefficients are expressed 
through the moments m at p, and solving differential equations Dpu = with 
specific numerical coefficients found in previous stages. Ultimately, \P p(x) = 
typj(x) is represented on each interval Aj as a linear combination of the basis 
solutions of Dpu = 0. 

Next, for each fixed x 6 Aj we consider equalities (3.4) for different j3 as 
a system of equations for the unknowns <fi (x),0 ; (x), with the known by 
now right hand side ^pj(x). This system of equations is a special case of 
the Prony system, as described in the previous section. Here the amplitudes 
aj are known to be ±1. (3.4) does not give the first equation in the Prony 
system, but it is just the sum of the amplitudes a,j, and we know it to be zero. 
So finally, we apply Procedure 2 and solve system (3.4), reconstructing the 
functions (x), 4>j t i{x) for each i 6 Aj. Now the functions <j) (x), (f>j t i(x) for 
j — 0, . . . , ft, I — 1, . . . , Sj, completely determine the domain G. Theorem 
3.1 is proved. 

As it was mentioned above, for the case of algebraic boundary segments 
explicit bounds can be given on the number of the moments required for re- 
construction. We provide these bounds in the case of polynomial boundaries, 
where the expressions are relatively simple and pretty sharp. 

Theorem 3.2 Under the assumptions of Theorem 3.1 let us assume addi- 
tionally that each boundary segment Sj is a graph of a polynomial y = ifij(x) 
of degree at most d. Then the domain G can be uniquely reconstructed from 
a collection of double moments 




(3.5) 



{m a> p : < p < 2(« - 1), < a < M(N, k, /3)} 



where 



M(d, k, (5) = max <^ 2{pN + 1)k - 2, (re + i)(pN + 1) 
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Proof: The functions ^ constructed in the proof of Theorem 3.1 are now 
piecewise-polynomials of degree at most /3N, and the conclusion follows from 
the corresponding result of [5] given in Section 2 above. 



4 Examples of Explicit Reconstruction 

The main object considered in this section is a compact plane domain bounded 
by a part of an elliptic curve 

y 2 = ax 3 + bx 2 + cx + d = f(x). (4.1) 

We assume that the roots X\ < x%, and x^ of the equation f(x) = are 
real, and that the function f(x) is positive on the interval (xi,X2). Hence, 
the equation y 2 = ax 3 + bx 2 + cx + d defines a compact domain Gel 2 . 

Consider a finite collection of corresponding moments 



Tn a> @ = / x a ■ y 13 dx dy. 
Jg 



Since the domain G is symmetric with respect to the axis Ox, we have 
m a,p — for odd values of /3, and 

2 f X2 

maw = -T, / x a ■ (ax 3 + bx 2 + cx + df +l/2 dx. (4.2) 

In order to simplify the formulae below, we introduce one more notation: 

and we call it "Moment". So, the equation (4.2) implies that 

M a ^ +2 = a ■ M a+3 , 2/ 3 + b ■ M a+2 ^ + c ■ M a+1 ^ + d ■ M a ^. (4.3) 

Our task is to determine the curve (4.1) from a finite (possibly minimal) 
collection of the Moments M a 2 /3- 

Let us calculate the Moments M a ^ for small values of the indices: 

M ,2 = a ■ M 3j0 + b ■ M 2j0 + c • M ifi + d ■ M , ; (4.4) 
M 1)2 = a ■ M 4j0 + b ■ M 3j0 + c • M 2 , + d ■ M 1)0 . 
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Hence, one can obtain a system of relations with unknown coefficients a, b, c, d, 
and to verify compatibility of the "data" {M Q ( g}. 

Here are two more methods of construction of similar relations: 

a) . Consider the Moment 

M Q)2 = / x a ■ (ax 3 + bx 2 + cx + df /2 dx. 

Since f(x\) = f(x 2 ) = 0, integrating by parts: 

du = x a dx, v = (ax 3 + bx 2 + cx + d) 3 ^ 2 , shows that: 

M a , 2 = - 2(q 3 +1) • (3a • M Q+3i0 + 2b ■ M a+2fi + c • M a+1 , ). 

In the cases a — 1, and a = we get 

• M 1>2 = 3a • M 4i0 + 26 • M 3i0 + c • M 2>0 . (4.5) 

and 

- o • M o,2 = 3a • M 3i0 + 26 • M 2 , + c • M lfi , (4.6) 
o 

b) . Similarly, since f(x\) = f(x 2 ) = 0, a simple change of the variables 
implies 

/ (3ax 2 + 2bx + c) ■ (ax 3 + 2bx 2 + cx + d)dx = 0, 

J X\ 

and hence, 

= 3a • M 2fi + 2b ■ M lj0 + c • M 0fl . (4.7) 



Consider now the system of three linear equations (4.5), (4.6), (4.7) with 
respect to the unknowns a,b,c. Let L 2 j[x 1 ,x 2 ] be the Hilbert space com- 
posed by corresponding functions defined on the segment [xi,x 2 ], endowed 
with scalar product: 



(F(x),H(x)} := / F(x) ■ H(x)y/f(xjdx. 



We have in these notations: M 4i0 = (x 2 ,x 2 ), M 3i0 = (x 2 ,^ 1 ), M lj0 = (x 1 ,^ ), 
M 2i0 = (x 2 ,x°) = (x^x 1 ), and M , = ' 
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So, the determinant of the system (4.5), (4.6), (4.7) equals 

/M 4 , M 3i0 M 2 , \ / (x 2 ,x 2 ) (x 2 ,x l ) (x 2 ,x°)\ 

6- M 3i0 M 2>0 M lj0 =6- (x 2 ,^) (x 1 ,^ 1 ) (x 1 ,^ ) , 
\M 2 ,o M 1>0 M ,o / V^ 2 ^ ) ( x ^ x °) ( x °> x °) J 

and this coincides (up to the factor 6) with the determinant of the Gram 
matrix of system of three polynomial functions x 2 , x l , and x°, which are 
linearly independent in the space L 2 j[xi, x?\. It is well-known that this 
determinant is strictly positive, thus, the system of linear equations (4.6), 
(4.7), (4.8) has a unique solution a,b,c. Then the coefficient d is uniquely 
determined from the equation (4.4), since we assume that M 0j o > 0. 
So, we have proved the following: 

Theorem 4.1. ([8]) In order to reconstruct an elliptic curve (4.1) it is 
sufficient to know 7 moments mo.o, TTti.O; rr^o, rri3,o, m^o, ^o,2 ^i,2- 

Note that such an "overdeterminancy" allows to obtain corresponding 
(nonlinear) relations between the moments listed above. 

Similar calculations illustrate theorems 3.1, 3.2 in a very simple case: 
{K, = 1, N = 2), see [7]. Here we reconstructed a triangle T C M. 2 with 
non-vertical edges from a given set of moments 

m ,o, mifl, m 2) o, m 3j0 , m ,i, my. 

Cf. also [20]. 



5 Invisible sets and functions 

Let V n denote the space of polynomials P(x\, . . . , x n ) and let a collection 
S C V n be fixed. We call a function / on M n iS-invisible if J Rn P(x)f(x)dx = 
for each P E S <zV n - A domain G C 1" is S- invisible together with its char- 
acteristic function. With obvious modifications this definition is extended to 
subsets of higher codimension and to distributions. 

In this section we discuss some examples of invisible sets and functions, 
coming from different fields. Besides Propositions 5.1 and 5.4 below, we 
do not provide new results, but rather an initial attempt to find similarity 
between several apparently unrelated problems. The common property we 
would like to stress is a remarkable "rigidity" and symmetry of the invisible 
objects. 
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5.1 S <zV n annihilating a fixed differential operator 

For a fixed partial differential operator T> in n variables it is natural to 
consider S C V n consisting of all P £ S with T>P = 0. 

5.1.1 Null quadrature domains 

Put T> = A to be the Laplacian, and denote Sh the corresponding set of 
harmonic polynomials Sh = {P £ V n , AP = 0}. 

A domain G C R™ is called a null quadrature domain if j G hdx = 
for all harmonic and integrable functions h. Taking h = P £ Sh we get 
a closely related notion, so null quadrature domains are essentially all the 
Sh- invisible sets. This class of domains has been intensively studied, and it 
has wide applications, in particular, in the investigation of the Newtonian 
potential, and of the filtration flow of incompressible fluid (see [2H, [251 EE] 
and references therein). 

Null quadrature domains include half-spaces, exterior of ellipsoids, exte- 
rior of strips, exterior of elliptic paraboloids and cylinders over domains of 
these types. It is known that in M 2 any null quadrature domain belongs to one 
of the categories above ([M]). A complete description of all null quadrature 
domains in higher dimensions has remained an open problem. A significant 
progress has been recently achieved (see [2H [25] and references therein) . 

5.1.2 Sets invisible for solutions of the wave equation 

Here we consider a somewhat artificial example which however illustrates the 
situation for another type of the operator T>. We put n = 2 and consider 
V = W = g^ty. in this case &w consists of all the polynomials P of the 
form P(x,y) = Q(x) + R(y). A function f(x,y) is SV-mvisible if and only 
if L a f(x, y)(Q(x) + R(y)) = for any polynomials Q(x), R(y). In turn, this 
is equivalent to the vanishing of all the moments 

f m;2 x k f(x,y)dxdy = f x k dx f f(x,y)dy and 

/ R a I l l f( x > V) dxdy = J y l dy J f(x, y)dx. 
This is equivalent to the identical vanishing of the functions 

F(x) = J f(x,y)dy = and H(y) = J f(x,y)dx = 0. So we have the 
following result: 

Proposition 5.1 f(x,y) is Sw -invisible if and only if F(x) = f f(x,y)dy — 
for each x and H(y) = J f(x,y)dx = for each y. In particular, this is 
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true for functions given by finite or infinite sums of the products (p(x)i/j(y) 
with f <p(x)dx = 0, f tj)(y)dy = 0. 

5.1.3 Vanishing conjectures of W. Zhao 

In a series of recent papers ([38, 39J and references therein) W. Zhao has stud- 
ied a number of vanishing conjectures which relate polynomials annihilating 
certain differential operators, invisible sets, and the well known Jacobian 
conjecture ([I]). 

For a given D, specifically, for T> = A being the Laplacian, the poly- 
nomials P have been considered satisfying the following condition: T> l P l = 
0, 1 = 1,.... This condition turned out to be closely related to the classical 
and generalized orthogonal polynomials. The following conjecture has been 
shown in [3S] to be equivalent to the Jacobian conjecture: 

Conjecture A If for a homogeneous polynomial P of degree four A l P l = 

0, 1 = 1,2,..., then A l P l+1 = 0, I > 1. 

It was shown in (38] that the vanishing of D l P l is equivalent to P being 
Hessian nilpotent — i.e. the Hessian matrix H(P) = ( 9 ®.q x . ) being nilpotent. 
In [39] Conjecture A has been closely related to the following Conjecture B: 

Conjecture B For a compact domain G C W 1 , for a positive measure \i 
on G, and for a polynomial P if all the moments J G P k dfi vanish, then for 
any polynomial q the moments J P k q dfi vanish for k ^> 1. 

In our language this conjecture can be reformulated as follows: if (G, y) 
is invisible for all the powers of P then it is "eventually invisible" for the 
sequence of polynomials P k q. Below we discuss this conjecture in somewhat 
more detail. 

5.2 Sets invisible by powers of a fixed polynomial 

In this section we discuss the vanishing problem for the moments J G P k d/j,, 

1. e. the conditions of invisibility of {G, n) for all the powers of P. Besides its 
appearance in Zhao's study of the Jacobian conjecture as above, this question 
is related to a wide spectrum of problems in Analysis, Algebra, Differential 
Equations, and Signal Processing. We shortly mention below only a very few 
of these remarkable connections. 
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5.2.1 One-dimensional case 

In one dimension the question is to describe all the univariate polynomials 
(Laurent polynomials, etc) P(x) and q(x) for which 



This question appears as a key step in understanding the classical Center- 
Focus problem of the Qualitative Theory of ODE's in the case of Abel equa- 
tion (see |10[ ITT] and references therein) . 

Even in this simplest case the answer (only recently obtained in [2?1 [29]) 
is far from being straightforward. In particular, it involves subtle properties 
of the polynomial composition algebra. To state the result we need the 
"composition condition" (CC) defined initially in [3\ and further investigated 
in pJJ [HJ [121 123 Ml [291 EQ] and in many other publications. 

Definition 5.1 Differentiable functions f(x) and g(x) on [a, b] C R are said 
to satisfy a composition condition (CC) on [a, b] if there exists differentiable 
W(x) defined on [a.b] with W(a) = W(b), and two differentiable functions 
F and G such that F(x) = ff f(x)dx and G(x) = J* g(x)dx satisfy 



If /, g are polynomials and they satisfy (CC) then W is necessarily also a 
polynomial. 

Composition condition implies vanishing of all the moments m,k (change 
of variables). Necessary and sufficient condition for vanishing of m,k for p, q 
polynomials is given by the following theorem: 

Theorem 5.1 ([271 EE]) The moments mk in (5.1) vanish for k = 0,1,... 
(i.e. [a, b] is invisible for P k q) if and only if q(x) = qi(x) H — • + qi(x), with 
I = 1,2 or 3, where qi, . . . , q\ satisfy composition condition (CC) with P(x) 
on [a, b], possibly with different right factors W±, . . . , Wj. 

Analysis of the case of rational functions, and, in particular, of Laurent 
polynomials (directly related to the Poincare Center-Focus problem for plane 
polynomial vector- fields) turns out to be significantly more difficult (see [28]). 




(5.1) 



F(x) = F{W{x)), G{x) = G{W{x)), x e [a,b]. 



(5.2) 
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If we allow P, q above to be only piecewise-polynomial (piecewise-rational) 
then another form of composition condition becomes relevant: a "tree com- 
position condition" (TCC) where W maps [a, b] not into R, but into a certain 
topological tree. Still under some restrictions a result similar to Theorem 5.1 
remains valid (see [12] ) . We hope that these recent developments can provide 
a better understanding of invisible D-finite domains, as above. 

5.2.2 Some examples in higher dimensions 

We start with a definition of a multidimensional composition condition (MCC) 
given in [16J, which directly generalizes Definition 5.1. (MCC) provides a nat- 
ural sufficient condition for the moments vanishing. However, as we shall see 
below, in n > 1 variables this condition is much stronger than the vanishing 
of the "one-sided" moments nik = f u F k (x)g(x)dx, k = 0, 1, ... . In fact, it 
is exactly relevant to the vanishing of the n-fold moments 



for all the nonnegative multi- indices a = («i, . . . , a n ). 

Let Q be an open relatively compact domain of M. n with a smooth bound- 
ary dQ. First we need for maps W : Q — > M n a definition generalizing to 
higher dimensions the requirement W (a) = W(b) in dimension one. 

Definition 5.2 ([IE]) A continuous mapping W : Q — > M. n is said to "flatten 
the boundary" dQ of Q if the topological index of W\sn is zero with respect 
to each point w eM. n \ W(dVt). 

Informally, W flattens the boundary dVL of Q if W(dQ) "does not have inte- 
rior" in M n . In particular, this is true if W\gn can be factorized through a 
contractible (n — l)-dimensional space X. The simplest example is when X 
is a point, so W mapping dfl to a point always flattens the boundary. We 
have the following simple fact: 

Proposition 5.2 (|16j) A mapping W : Q — > W 1 flattens the boundary dQ if 
and only if the integral J n H(W(x))dW(x) vanishes for any function H{W). 

Now let Fx , . . . , F s be differentiate functions on Q and let fi be a measure 
on Q given by its density g(x): dfi(x) = g(x)dx. 




(5.3) 
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Definition 5.3 ([IS]) Functions J 7 }, I = l,...,s and a measure /i on fi 
satisfy multi- dimensional composition condition (MCC) if there exists a dif- 
ferentiate mapping W : fi — > ~R n , flattening the boundary <9fi, functions 
Fi(w), I = 1, . . . , s, and g{w) on W 1 such that Fi(x) = Fi(W(x)), I = 1, . . . , s, 
and d/j(x) = g(x)dx = g(W(x))dW. 

The following simple proposition (implied directly by Proposition 5.1) shows 
that (MCC) is sufficient for moment vanishing: 

Proposition 5.3 If a function F and a measure \x on fi satisfy (MCC), then 
all the moments rrik = J n F k (x)g(x)dx, k = 0, 1, . . . , vanish. 

Consider the following example: let fi C M™ be defined by P(x) < 1 for a 
certain polynomial P( ) E W 1 . For each j = 1, . . . , n define 

Sj to be a collection of polynomials Sj = {Q n (P)§^-} with Q an arbitrary 
univariate polynomial. 

Proposition 5.4 For each j — l,...,n, fi, Q(P), d/j = §^-dx satisfy 
(MCC), so the domain Q is invisible for Sj. 

Proof: Define W : Q -)■ R n by W(xi, . . . ,x n ) = (yx, . . . , y n ), with y { = 
x ii i 7^ i) Vj — F(xi, . . . , x n ). W maps the boundary dVL into the hyperplane 
{ Vj = 1} c R n , so it flattens dtt. Now, Q(P) = Q{W), where Q(y u ...,y n ) = 
Q(Vj), and 

dP 

d/j, = — — dx = dxi ■ ■ ■ dxj_i ■ dP ■ dxj + \ ■ ■ ■ dx n = dW. 

OXj 

Let us now describe a situation where (MCC) is a necessary and sufficient 
condition for invisibility. Consider double moments of the form 

m k ,i = I P\x, y)Q\x, y)r(x, y)dxdy, k, I = 0, 1, . . . , flcl 2 . (5.4) 
Jn 

We shall assume that P in the domain of consideration satisfies 7^ and 
consider Q of the form a < P(x,y) < b, c < y < d. The functions P,Q,r 
in (5.4) are assumed to be real analytic, and Q is assumed to have a simple 
critical value on each level curve of P inside Q. 

Theorem 5.2 Q16J) Under the above assumptions all the moments mkj, k, I = 
0,1,..., vanish if and only if Q, P, Q, rdx satisfy (MCC) . 

So the domain Q as above is S'-invisible for S consisting of all the products 
P k Q l r if and only if fi, P, Q, rdx satisfy (MCC). 
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5.2.3 Mathieu conjecture and Laurent polynomials 

Zhao's Conjecture B above has been motivated, in particular, by the following 
conjecture of O. Mathieu Q26J), closely related to many important questions 
in Representation Theory: let M be a compact Lie group. Denote F(M) the 
set of M-fmite functions on M (i.e. polynomials in all the characters on M) 
and let fi be the Haar measure on M. 

Conjecture C. If for some f(x) 6 F(M) 

/ f k {x)dfi{x) = 0, k = 1,2, ... (5.5) 

then for any g(x) G F(M) we have J M f k (x)g(x)d/j,(x) = 0, k 3> 1. 

This conjecture is known to imply the Jacobian conjecture Q26J). In our 
language, it states that if M is invisible for f k , it is eventually invisible for 
f k g with any M-fmite function g. 

Conjecture C has been verified in [T5] for the Abelian M, i.e. for M 
being the n-dimensional torus T n . In this case M- finite functions are Laurent 
polynomials in z = (z±, . . . , z n ), Z{ G C, \zi\ — 1. In fact, the following result 
has been established in [T5] : 

Theorem 5.3 Let f(z\, . . . , z n ) be a Laurent polynomial. Then the constant 
term of f k vanishes for k = 1, 2, . . . if and only if the convex hull of the 
support of f does not contain zero. 

Here the support of / is the set of multi-indices of all the monomials in 
/ with nonzero coefficients. Theorem 5.3 immediately implies Conjecture 
C since under its conditions the support of f k eventually gets out of any 
compact set on Z n , in particular, out of the support of g. 

Recently a rather accurate description of moment vanishing conditions 
for one-dimensional rational functions and, specifically, for Laurent polyno- 
mials has been obtained in [2H]- In particular, an extension of the result of 
Duistermaat and van der Kallen ([IS], Theorem 2.1 above) obtained in [2B] 
provides such conditions: 

Theorem 5.4 ([28], Theorem 6.1). Let L(z) and m(z) be Laurent polyno- 
mials such that the coefficient of the term - in m(z) is distinct from zero. 
Assume that J sl L k (z)m(z)dz = 0, k ^> 1. Then L(z) is either a polynomial 
with zero constant term in z, or a polynomial with zero constant term in - . 
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As it was explained above, this property implies that j sl L k (. 
0, k ^> 1 for any Laurent polynomial h(z). In particular, we get 



(z)h(z)dz 



I 



L k (z)g(z)m(z)dz = 0, k > 1 



for any Laurent polynomial g(z). Therefore Zhao's Conjecture A holds for 
S 1 and the measure d/i(z) = m(z)dz. 

In [30] under a stronger assumption of vanishing of the moments starting 
from the initial indices, we get the same conclusion assuming that only a 
"horizontal strip" of the moments vanish. 

5.3 Complex moments 

Problems of reconstruction of sets and functions from their complex mo- 
ments, and, in particular, vanishing conditions for complex moments form 
an important field of investigation in Several Complex Variables, in Inverse 
Problems in PDE's and in related fields. We give here only a few examples 
illustrating connections with our setting. 

5.3.1 Wermer's theorem and later developments 

The classical theorem of Wermer ([3?]) gives conditions for vanishing of all 
complex moments j^x l y^dx for a closed curve 7 C C 2 : this happens if and 
only if 7 bounds a compact complex one-chain. See [131 1221 1231 EZ| and 
references therein for an accurate statement and further developments. In 
our terms 7 is invisible for all complex moments if and only if it bounds a 
compact complex one-chain. 

The theorem of Dolbeault-Henkin ([13]) gives a remarkable extension of 
Wermer's condition to the case of curves 7 bounding a compact complex one- 
chain in the projective space. In particular, in such case the moment generat- 
ing function satisfies a non-linear Burgers-type partial differential equation. 
This last fact can be reinterpreted as an invisibility of 7 for certain combi- 
nations of the complex moments. 

Let's assume now that 7 C C 2 is an image of a not necessarily closed 
curve a C C under a rational mapping (P, Q) : a — > C 2 . In this case a more 
accurate form of Wermer's theorem can be obtained: 
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Theorem 5.5 (|28j. Theorem 5.2) The moments m^j = j^x l y^dx vanish 

for j,j > 1 if and only if there exist rational functions P,Q,W such that 
P(z) = P(W(z)),Q(z) = Q(W(z)), the curve a = W(a) C C is closed, and 
all the poles of P, Q lie on one side of the curve a. 

In particular, if the moments rriij vanish for i,j 3> 1, they in fact vanish 
for all i,j > 0. 

In our terms, if the curve 7 = (P, Q)(o~) is eventually invisible for the complex 
moments m»j then it is closed, it is an image under P, Q of the closed curve 
a C C, and 7 bounds the compact complex one-chain (P, Q){G) where G C C 
is the domain bounded by a and free of poles of P, Q. 

5.3.2 Complex moments of planar domains 

This problem has been intensively studied, in particular, in relation with 
the filtration flow of incompressible fluid, and with the inverse problem of 
two-dimensional Potential Theory (see I2H ES] and references therein). 
In particular, in [36] one can find a discussion of the non-uniqueness of re- 
construction. 

For reconstruction of polygonal domains see [171 El]- m general, an impor- 
tant class of quadrature domains (see [21 [21] and references therein) provides 
a natural framework for a study of the "finite dimensional" reconstruction 
problem, as well as its possible non-uniqueness, in particular, the invisibility 
phenomenon. 

6 Conclusions 

The results of the present paper leave open some important questions: 

1. We have insisted on a requirement that the singular points of the 
differential operator D are outside the projection of the domain to be re- 
constructed. This requirement excludes some important classes of domains. 
In particular, it prevents reconstruction of general semi-algebraic domains G 
on the plane. Indeed, typically the projection of G to the x-axis will have 
singularities on the boundary dG. These singularities will be necessarily also 
singularities of the differential operator D annihilating the corresponding 
algebraic functions. The construction of [5] allows for an adaptation to sin- 
gular situations: just, on the last step we have to use the bases of "singular 
bounded solutions" of D = 0. See, for example, [T]. 
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2. An important parameter entering the procedure of reconstruction of a 
D-finite function / is the maximal number /x of its initial moments that can 
vanish, unless all the moments vanish identically (see an example in Section 
2.1 above). Recently we've shown that this number depends only on the 
combinatorial data in case where singular points of D differ from the jump- 
points of /. If some of the singularities of D are the jump-points, we expect 
an explicit bound on \i through the size of the coefficients of D. 

3. Robustness estimates. A "quantitative version" of the question of 
bounding of the number /i is to bound all the moments of / through its 
first fi moments ("moments domination"). We expect this problem to be 
a central one for the robustness estimates of the reconstruction procedure. 
Indeed, moments domination implies directly a bound on the norm of / 
itself through its first /i moments. This question is also directly related to 
the analysis of periodic solutions of Abel equation ([TT]). 
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